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Abstract: 

We discuss a possible generation of color suppressed B-decays amplitudes through 
a soft final state interaction. As a typical example, we consider in detail the decay 
B° — > D°7r° (and also B° — > 27r°). We show that in the approximation of the two 
particle unitarity and at zero order in a s this process can be related to the weak 
decay B° — > D + n~ followed by the strong charge exchange scattering in the Regge 
kinematics. We estimate the amplitude of this process using the light cone QCD 
sum rule technique and find that it is supppressed as a power of 1/mg in comparison 
to the amplitude generated by the effective non-leptonic Hamiltonian, but remains 
important for the physical value of 
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1 Introduction 



Understanding final state interactions (FSI) effects in weak decays is a challenging 
task of acute interest, both theoretical and phenomenological. Among plenty of 
problems, strongly influenced by possible FSI, one can mention a determination 
of the Kabayashi-Maskawa matrix elements, weak CP violating phases, etc. It is 
tremendously difficult to address the issue of FSI on a quantitative level. Still, one 
may believe that with going to heavier (c,b) quarks possible FSI would become less 
and less important. In most intuitively appealing form this expectation has been 
expressed by Bjorken [3J, who viewed this behavior as a typical manifestation of 
the color transparency phenomenon This scenario can be explained as follows. 
When a virtual W-boson decays into a ud pair with a low invariant mass, this 
pair is initiallly only weakly interacting with the bulk of the system as its color 
dipole moment is small. On the other hand, the pair escapes far away from the rest 
of decay products before starts to interact strongly, as a result of the relativistic 
time dilation. It follows on these grounds that in decays like B — > Dtt FSI die off 
asymptotically as mj, — > oo. In effect, weak non-leptonic decay amplitudes factorize 
into a product of hadron matrix elements; hard gluons, which otherwise would 
violate this factorization, do not manage to appear on time. 

It is worth noticing that this scenario, while is very attractive in view of its sim- 
plicity and appeals to the one's intuition, is only qualitative and does not predict the 
rate of vanishing FSI corrections when — > oo. Moreover, by its very appearence, 
this argumentation is likely to be applicable only to FSI proceeding via a hard glu- 
ons exchange. In fact, the factorization property has been proved to all orders of 
the perturbation theory in the heavy quark limit || m c ,m;, — > oo , m c jm\, fixed. 
On the other hand, in the case of non-perturbative soft FSI the time dilation is not 
so important, and one may expect that those effects could be sizeable in spite of the 
large b-quark mass. Moreover, if initial wave packages hadronize without gluon FSI 
when their separation is not yet large ( ~ 0.3 fm ), then there usual hadron-hadron 
scattering would be a dominant mechanism of FSI. If this kind of FSI is strong 
enough, it presumably could screen the effect of absence of perturbative FSI, which 
is due to the color transparency phenomenon. The present paper is just an attempt 
to provide a semi-quantitative estimate of a possible role of hadronic FSI corrections 
in non-leptonic B-decays of the type B — > Dn. Our approach is essentially based 
on the existence of large parameter in the problem ( mg ) , and differs from more 
traditional approaches to FSI in weak decays. 

We would like to start our presentation by explaining why FSI corrections are so 
elusive and hard to treat in most non-perturbative methods which have been applied 
to non-leptonic B-decays. As has been first explained in Ref.Q, the problem lies 
in the fact that physical amplitudes are initially extrapolated into the Euclidean 
region where a calculation is feasible, and then are analytically continued back into 
the Minkowsky space-time (this is the typical strategy e.g. in the QCD sum rule 
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method 0] or lattice calculations). Atter the Wick rotation fastly oscillating pole 
terms (corresponding to particles interacting in the final state) become falling ones 
and die off in the deep Euclidean region. On the other hand, when the amplitude is 
restored from a dispersion relation, one integrates actually over some finite region 
of the invariant mass where these oscillating terms are averaged to some smooth 
function. In effect, amplitudes are always real within these methods; scattering 
phases cannot be calculated. An analogous statement in the lattice approach has 
later become known as the Maiani-Testa no-go theorem |6|, [7]]. 

In this paper we try to estimate semi-quantitatively FSI corrections to the B- 
decays of the form B — > Dir. We consider the limit of the heavy b-quark. As 
for the c-quark, we distinguish between two different limits: (i) c-quark is heavy 
with nib ^ m c and (ii) c-quark is light, the corresponding process is B — > tttt. To 
simplify the problem and avoid complications due to presence of many channels, we 
have chosen to look for FSI effects on color suppressed B-decays (more specifically, 
we consider a particular case of the B° — ► D n decay) and work at zero order 
in a s . Our approach is strongly correlated with the fact that the energy excess in 
B-decays of this type is large which, as we show, imposes a relevance of the Regge 
approach. Surprisingly, we find that though FSI effects are suppressed as a power 
of 1/me in comparison with direct amplitudes, they remain to be of the same order 
of magnitude at the physical value tub = 5.27 GeV. Here we have to mention 
a recent paper || where similar to ours ideas have been applied to another class 
of B-decays. We still believe that we manage to reach more concrete motivations 
and results for our particular decay mode. We will compare our conclusions to 
those of Ref. || in the last section of this paper. In Sect. 2 we show that in the 
approximation of two particle unitarity the decay of interest can be related to a 
process of strong charge exchange scattering in the Regge kinematics and calculate 
a discontinuity of a corresponding Feynman diagram. In Sect. 3 and 4 we apply the 
light cone QCD sum rule technique (see e.g. |L(J for a mini-review) to estimate 
parameters of interest. As a by-product of our study, we derive new light cone 
sum rules for the B — ► D transition form factor (Sect. 3) and pD coupling constant 
9 p dd (Sect. 4) . In Sect. 5 we compare the calculated reggeon-induced amplitude 
to that proceeding directly via corresponding terms in the effective non-leptonic 
Hamiltonian. Our result is consistent with the color trasparency picture as long as 
the former is suppressed approximately as 1/mg in respect to the latter. However, 
we find that FSI corrections are pretty large for the physical value tub — 5.27 GeV 
and constitute 40 - 100 % of the amplitude of direct process. Sect. 6 contains a 
discussion of our results. 
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2 Discontinuity of the B —> D°tt amplitude 



We consider the following simplified model for the process of interest. Since at 
the zero order in a s the direct decay B — > D 7t° does not occur, we seek for a 
hadronic intermediate state which would proceed non-perturbatively into the given 
final state. The most obvious candidate is the D + n~ intermediate state which is 
related to D°7r° by the quark rearrangement. In what follows we restrict ourselves 
by this particular contribution only, i.e. work in the approximation of two particle 
unitarity relation. Corrections due to multiparticle intermediate states of the form 
D + + rm are presumably small, though accurate estimates are troublesome. Still, 
calculations within the quark gluon string model 1T1J indicate that the averaged pion 



multiplicity is low (2-3), and thus we hope that the two particle approximation is 
satisfactory with an accuracy not worse than the factor two. 

The discontinuity of the amplitude is given by the expression 



DiscM(k 2 ) = /^^(p-A;) 2 )^ 2 ^)^ 2 ,^-/-^ 2 )^^^^-), (1) 

where T(k 2 , {k — I — p) 2 ) is the strong charge exchange amplitude for the process 
D + (p) + Tr~(k - p) -> D°(k - I) + 7r°(Z) and A W (B -> Dtt) is the weak amplitude 
for the transition B° — > D + {p) + 7i~(k — p). Using the bare weak Hamiltonian 

H w = ^ Vcb V: d (cT,b) (dT.u) (2) 

( r M = 7^(1 — 75) ) and taking into account delta functions in the above formula, 
we obtain 

A w (B(k) - D(p)n(k - p)) = if^ Vcb V: d (k 2 - m 2 D )f + (0 7 k 2 ) , (3) 

where the form factor /+(0, k 2 ) is defined according to 

(D(p)\cT,b\B(k)} = f + ((k-p) 2 ,k 2 )(k,+p,) + f-((k-p) 2 , k 2 ){k,- Plx ) (4) 

Here we have introduced an additional variable k 2 into the standard definition of 
the B — ► D form factor bearing in mind a subsequent analytical continuation in k 2 . 
We note that on the mass shell the form factor f + (0,m 2 B ) can be expressed via the 
Isgur-Wise function [J2j 

f (n ^ m B + m D (m B -m D ) 2 

/+( °' mB) = + 2m B m D } (5) 

We will not use, however, this correspondence but rather derive a new light cone 
sum rule which is more suitable for an analytical continuation and refers directly 
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to /+(0, m B ) . We delay a corresponding calculation until the next section and 
concentrate first on the structure of the integral (1). An integration of Eq.(l) yields 

Di8cM(s) = ^U^=V cb V: d f + (0,s) £ d(cos 9)T(s, t) {s ~£ d) * Q(s - mg,) (6) 

Here 9 is an angle between the momenta k — I and p : 

t = {k-l-pf = ~ {k2 ~™ l) \ l ~ cos 9) (7) 

The following observation is crucial for the rest of this paper. First, we note that 
Eq.(6) contains the s-wave partial strong amplitude, in accord with the angular 
momentum conservation. Therefore, one could think at the first sight that the Regge 
phenomenology, as being stemming from a resummation of high partial amplitudes, 
is of no use here. However, it has been known for a long time that partial amplitudes, 
obtained from Regge amplitudes by extrapolation into a low energy region, reveal 
resonance structures. This means that a correct answer is to be obtained by first 
resummation of partial waves into a Regge form, and then re-expanding it back into 
partial amplitudes. It is useful for a moment to think of the final state interaction 
as being proceeding through a resonances formation. Anticipating the later result 
(see Eq.(13) below) we note that the discontinuity in the dispersion integral is a 
slow (power-like) function of s. Therefore in view of a large parameter m 2 B in the 
denominator a very large region of s from the threshold s = (mjj + m,) 2 to s ~ m 2 B 
contributes in the dispersion integral. It is quite clear that in this parametrically 
large region many overlapping resonances are created. According to the well known 
particle - reggeon duality phenomenon s-channel resonances can be be described by 
a t-channel reggeon exchange To avoid possible misunderstanding, we should 



emphasize a difference of our case from the duality picture which applies, strictly 
speaking, only to an imaginary part of scattering amplitude averaged over some 
region in energies. Our consideration applies directly to the scattering amplitude, 
as follows from the discussion above and Eq.(6). Thus, we conclude that the integral 
in Eq.(6) is dominated by a strong scattering amplitude in the Regge kinematics. 
Taking into account the leading p— meson trajectory, we obtain 

T(s,t) = P(tm(-) a(t) (8) 

Here a(t) — a(0) + a't is the Regge trajectory, s is a scale parameter, £(£) is the 
signature factor 

m = j-^H^*)] (9 ) 

and j3{t) stands for the reggeon residue. In what follows we adopt a model for the 



functional t-dependence in (3{t) obtained in dual models JT3[ and the quark-gluon 
string model |TJj : 

m = (io) 

T[a{t)\ 
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We next calculate the integral in Eq.(6) neglecting the t- dependence of the product 
T(a(t)) sin 7rcx(t) 



T(a(t))smna(t) ~ r(a(0)) sin7ra(0) 
and approximating t ~ — s/2(l — cos 0) . Then 



;n) 



d{cos6)T(s,t) 



2/3(0) 



«(0) 



1 — exp [— a's In — ] 1 — exp [— a's(ln — — wr)] 



In -2- 

«0 



In-*- 

«0 



(12) 



Since the exponents are small in comparison to unity above the threshold s = w 
we finally obtain 
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Disc M(s) 



^%V cb V: d Uf + {^ s)F(s)^ r (-V (0) &(s - ml) (13) 
Sir a/2 gnat's Vs / 



where 



F(s) 



1 



+ 



(14) 



In (2a 1 s) ' In (2a's) - m 

In these formulae the intercept a(0) ~ 0.5 and the slope of the Regge trajectory 
oi ~ 1 Gev~ 2 . Furthermore, it is important to stress that the scale parameter 
so = Sq D is not constant, but varies with m c . To parametrize this dependence, we 
will use the relation obtained in the quark gluon string model [fbl 



,7rD 



S -d — 



a: 



x 



D 



2A 



(15) 



where (xjj) is the average momentum fraction carried by the d-quark in the ir- 
(D-) meson and A ~ 300 — 400 MeV is the average virtuality of the light quark in 
the D-meson. For the physical value ~ 1.9 GeV one obtains Sq D ~ 2(a')~ 1 . It 
is worth noticing that the standard dual model arguments would rather suggest a 
universal expression for the scale factor Sq D ~ (a/) -1 . 



3 Light cone sum rule for the B —> D form factor 



Before proceeding to a calculation of the dispersion integral we need an estimate for 
a functional s dependence of the B — > D transition form factor /+(0, s) . The idea 
of a method presented in this section is borrowed from a very similar problem of the 
pion form factor at intermediate momentum transfers, which is most conveniently 
solved by the light cone QCD sum rules method ]T3| (see also Ref. [0]). We consider 
the following correlation function 
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T^fak) = % J dxe ikx (D{p)\T{cj^b{x)hul 5 d{0)} |0> (16) 

for Minkowskian k 2 = and large Euclidean momentum (p + k) 2 — > — oo. The 
B-meson contribution to this correlation function is 

- w ' ^W /t(t ' i+ - • < i7) 

where ellipses stand for other Lorentz structures. On the other hand, the correlation 
function (12) can be calculated in QCD at large Euclidean momenta (p+k) 2 . Indeed, 
in this case the virtuality of the heavy quark, which is of order ml — (p + k) 2 , is large. 
Thus, we can expand the heavy quark propagator in powers of slowly varying fields 
residing in the D-meson, which act as external fields on the propagating heavy quark. 
The expansion in powers of an external field is also the expansion of the propagator 
in powers of a deviation from the light cone x 2 ~ 0. The leading contribution is 
obtained by using the free heavy quark propagator in the correlation function (12). 
We obtain 

T^ip, k) = J (27r)4K _ /2) e^(D(p)|c(x) 7 ,(/^ + m b ) lulb d(m (18) 

In this formula a path-ordered gauge factor between the quark fields is implied, as 
required by gauge invariance. In the particular case of the Fock-Schwinger gauge 
x^A^ix) = this factor is equal to unity. We see that the answer is expressed via the 
one-meson matrix element of the gauge invariant non-local operator with a light-like 
separation x 2 ~ 0. This matrix element defines a light cone D-meson distribution 
amplitude. The leading twist distribution amplitude <f>(u) is introduced as follows 

(D(p)|c(x)7^7 5 d(0)|0) = -ifoPn / du exp (iupx)(f)(u) (19) 



Within this definition we obtain to the leading twist accuracy (hereafter u — l — u) 

T^(p, k) = -if D Pnh f du —J- — _ / — "Mi + " ' ( 20 ) 

jo + uum D — u(p + k) z 

To match Eq.(19) with the B-meson contribution (16) to the correlation function 
(15), we note that Eq.(19) can be re- written as the dispersion integral with (m 2 + 
uum 2 D )/u being the mass of an intermediate state. The duality prescription tells 
that this invariant mass has to be restricted from above by the duality threshold 
so = (rrib+mo) 2 , where m does not depend on m&. For physical value m& ~ 4.7 GeV 
one gets so ~ 35 GeV 2 , thus m ~ 1.2 GeV. We then obtain the following sum rule 

2 2f , B , fr /+(0, m B ) = fn f 2mo du _ , - V2 +■■■ (21) 

m B — (p + qy m t + uum D — u(p + q) z 

To analyze a parametric behavior of the form factor / + (0,ms), we first note that as 
ml 3> ml we can neglect the second term in the denominator in the right hand side 
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of Eq.(20). Next one has to distinguish between two different kinematical regimes : 

(i) In the limit m& —>■ oo and m c —>■ D-meson becomes pion, the corresponding 
process being B — > 7r + 7r~ ' -> 2tt°. In this case the asymptotic as w — > I behavior of 
4>{u) is ~ 6(1 — u) and then 

/+(0, m B ) ^ -^2" ~ m B ( 22 ) 

(here we have used the fact that for mg — > oo /q ~ mg up to logarithmical 
corrections.) This dependence is well known from the light cone sum rules for the 
B — > 7r transition 0, [Tcfl . 

(ii) More interesting is the case m c , — > oo with mj, ^> m c . Then the D-meson 
distribution amplitude is strongly peaked at equal velocities of quarks, i.e. concen- 
trated near u — 1 [IT]]. We adopt the parametrization [ IB 



d (m) = A D — — — — , (23) 

where e = A/(A + m c ) and Ad is fixed by the normalization ( Ad — 0.67 for 
m c ~ 1.3 GeV and A ~ 0.33 GW ). Using Eqs. (21) and (23), we obtain 

/ + (0,m B ) - f ADfDm i . 1 ■ - ^. 5/2 (24) 



'lm, 



More faster in comparison with (22) fall off of the form factor / + (0,mg) indicated 
in Eq. (24) appears quite reasonable. As 4>d{u) is strongly peaked near u = 1, 
the integral in Eq. (21) becomes very sensitive to a change of the integration limit. 
Note in passing that this is precisely the reason why a naive ansatz of the type 
4>(u) = 8{u — 1 + A/m c ) would produce completely wrong answer. 



4 Residue factor (3(0) and g P DD coupling constant 



In this section we obtain an estimate for the residue /3(0). To this end, we make 
use of the fact that a reggeon-particle vertex is an analytical continuation of the 
t-channnel resonance coupling constant into the region t ~ |13] . We thus consider 



the charge exchange process D + n — > D°ti proceeding for t ~ m 2 p via the p- meson 
exchange. The asymptotic ( for s — > oo ) expression for the amplitude is 

2s 

T(s,t) = gpnnQpDD- ? (25) 

t — trip 

The p — 7r vertex is well known: 



— Tfl 

9pn7T ^ 9 P , (0|e u «'y M i4 + e d d^^d\p) = — ^ , (26) 



9 P 
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( Bp stands for the p- meson polarization vector) and thus we have to estimate only 
the pD vertex g P DD- Here we again use the light cone QCD sum rules technique 0. 
Consider the correlation function 



T^(p, q)=ijdx e^(p(g)|T{rf7 M 7 5C (x)c 7 ,75M(0)}|0) (27) 
The contribution of interest is 

+ <I,)(ep)9 pDD f D {m 2 D _^l-(p + q )2) ^ 

On the other hand, we calculate this correlation function at large Euclidean momenta 
p 2 , (p+q) 2 — > — oo using the Operator Product Expansion near the light cone P, [TO | . 
To the leading twist accuracy we obtain 

T =-pq (ep) f m f du jM^M ( 2 g) 

where the twist 2 distribution amplitudes g±\u) and &\\(u) are defined as follows 
T% HI: 



(p(q)\d{x)^u(0)\0) = q p {ex)f p m p f due iqx $\\(u) 

+ e p f p m p f 1 due tqx g { l\u) (30) 

J 



(p(g)|d(a;)7 M 75«(0)|0) = ~f p m p E puXa E u q x x a due iqx g { f> (u) (31) 
To the lowest conformal spin accuracy the distribution amplitudes have very simple 



forms [19 



3 

$H (u) = i-uu(2u - 1) 

g { l\ u ) = 6uu (32) 

We next use the identity (p + uq) 2 = up\ + up 2 — uuq 2 with p\ = p 2 and p\ = (p + q) 2 
and the formula for the double Borel transformation in respect to p\ and p\ 



BM ^ [ml + uuq 2 -up 2 -up 2 Y = ' *** { ~lt 2 ' WTM 2 ^ " U ° 

(33) 



where 

M 2 M 2 Mi 
Ml + Mf ' U °~Mj+M'i 

In a final expression we set M\ = M 2 = 2M 2 and account for the fact that $ii (1/2) 
0. We then arrive at a simple sum rule 



f p m p (a) 1 / m 2 c m 2 s \ 

QpDD - -^fsY^) (exp {- w - ^) - exp (-— ) j (35) 
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Here the second exponent corresponds to a continuum subtraction with s ~ 6 GeV 2 . 
This sum rule yields 

9 p dd ^ 4 (36) 
On the other hand, the Regge formula taken for t > provides 

nM ^ a «3>_^_ (37) 

Comparing these expressions and taking into accout Eq.(lO), we arrive at an estimate 

(3(0) ~ 2irT{a(m*)]g pn7T g pDD ~ 140 (38) 



5 Reggeon versus direct amplitude 



To obtain the amplitude of the process of interest, we write down an unsubtracted 
dispersion relation 

. « 1 f°° , . Disc M(s') , . 

Ms = — / ds' — — 39 

v 1 2mJ m % s'-s-ie v 1 

To calculate the integral, we use a standard trick. The integral is first calculated 
for negative s = —m 2 B and then one makes the analytical continuation m 2 B — > 
m 2 B exp (— in) to positive s = m 2 B . To make subsequent formulas feasible, we neglect 
in this calculation the s-dependence of the function F{m B ) (see Eq. (14)). It is 
convenient to parametrize this quantity in the form 



F(m B ) = \F\exp(i<f>) (40) 



with ( a = In (2a'm 2 B ) ) 



i 

r 2 \ 2 



a \ 7r 2 + a 2 y 
a 2 

= atan^ ~ 50° (41) 

a 2 + it 2 — na 

(here numerical values refer to the physical mass ~ 4.7 GeV). Again, we distin- 
guish between two different limiting cases: 

(i) mi, — ► oo and m c — > 0. The corresponding physical process is i? — ^ 7r + 7r~ — > 27T°. 
In this case the scale parameter is the standard s^ = (a') -1 . The asymptotic with 
m h — > oo behavior of the -B — > 7r transition form factor is given (in the chiral limit) 
by Eq. (22). However, the asympotic estimate (22) cannot be analytically continued 
down to the threshold value s = 0. To get an interpolating formula, we note that 
at mj — > /+(0, mg) reduces to the pion electomagnetic form factor which equals 
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unity at zero momentum transfer. We model this by the interpolating formula which 
is approximately valid for all in the dispersion integral 

U(0,s)= 6f ^f 2 (42) 

Neglecting for simplicity the m B -dependence of during the calculation of the 
integral, we finally obtain 

M£(m B ) = ^Wa^B ~ W (43) 
V2 8ttJti a' 3 ™b/b 



(here the subscript "r" means the reggeon-exchange induced mechanism of weak 
decay.) 

(ii) m c , rrib — > oo with rrib ^> m c . The corresponding physical process is B — ► 
D + n~ — > D 7i°. Using Eqs. (15), (24) and (38), we obtain in this case 

AffiH = %Vt^ elFlm^e^/^-W (l-^L) (44) 



7T 



In deriving this formula we have restricted ourselves by the linear order in n 



"^Dl m B- The answers we have found are to be compared with perturbative QCD 
amplitudes of the direct decays which are induced by the well known effective non- 
leptonic Hamiltonian (we display only the Cabibbo favored part) 

H ef f = ^WfcfrOOO? + CMO u 2 ] (45) 

where 

0\ = (cV^idV^u) and O u 2 = (Jl»(cl» (46) 

The Wilson coefficients Ci{fj) are due to the renormalization of the bare Hamiltonian 
Hw ~ 0\ by hard gluons with virtualities larger than fi 2 = 0(m|). In the leading- 



log approximation with fi ~ 5 GeV , nfi — 5 and A^jg ~ 200 MeV [gC 



Ci = 1.117 and C 2 = -0.266 (47) 
For the B° — ► D°7r° decay O2 factorizes and 



(DVl^/zlS ) ~ (C 2 + ^)(L>|cr^|0)(7r|dT M 6|S) 

+Ci( J D7r|2(cr M t a M )(Jr M t a 6)|5) (48) 

The non-factorizable contribution to this decay which is given by the second matrix 
element in this formula has been estimated in Ref. [J2T]] by the light cone QCD 
sum rule method. There has been found the following estimate for the ratio of the 
non-factorizable to the factorizable 1/N C amplitudes : 

r - M(nf) AT C Pa (rr\dG a ^ 5 b\B) ^ 

~ M(f) ~ A^p D {ml - ml)f+(ml) ~ ' 1 1 
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where the the B tt transition form factor /+ is defined as in Eq.(4). The direct 
decay amplitude is then 



M^ dir) i 



ni 



B) 



73 



fDf«{m 2 D ){m 2 B - m 2 D ) ~ m 



i /2 



(50) 



We therefore see that the reggeon- induced amplitude is suppressed as l/m% in re- 
spect to the direct one, while at m# = 5.27 GeV one obtains 



M 



(r) 



0.09 



C 2 + f (l + r) 



3 j(0+7r/4-2?7/7r) 



0.18 + 0.58i 



(51) 



The last equality in this formula is obtained with the particular value (49). Repeat- 
ing the same procedure yields a similar formula for the B — > 2tt° decay. In this 
case (if one neglects a penguin contribution) the amplitude factorizes analogously 
to (48), but this time the suppression is only 1/itlb- Numerically one finds in this 

case 

M ^ °- 12 B i(^/4) (52) 



Mi d : r) 



where the parameter r is defined analogously to Eq.(49). Our final results Eqs.(51) 
and (52) suggest that the rescattering amplitudes are not small and must be added 
to the amplitudes of direct decays. In this sense the latter can be thought of as 
radiative corrections to "tree" level amplitudes induced by soft FSI. We conclude 
that perturbative QCD most likely does not exhaust non-leptonic B-decays. The 
physical b-quark mass, albeit large, is still beyond the onset of the asymptotic 
QCD regime in nonleptonic heavy meson decays. An example of a similar interplay 
between soft and hard QCD is discussed in the next section. 



We would like to end up this section with a comment on a relation between the 
direct and rescattering amplitudes. For definiteness, we will focus on the QCD sum 
rule approach. As has been explained in the introduction, any information about the 
imaginary part of the amplitude is lost there. Our method therefore yields the entire 
phase of the decay. The issue of the real part of the rescattering amplitude is slightly 
more subtle. Real parts of the amplitudes calculated in this paper correspond to the 
so-called "parasite" contributions to the sum rules (cf. Ref. f|]) and are due to a 
possibility to produce states below the B-meson pole in the corresponding channel. 
If the latter are omitted from an analysis of QCD sum rules (this require a careful 
analysis of corresponding Operator Product Expansion), then the two amplitides 
are simply added. 



6 Concluding remarks 



In this paper we presented a semi-quantitative calculation of FSI effects in color 
suppressed B-decays. Our results suggest a rather dissapointing conclusion that the 
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soft physics prevents currently used methods from reaching an accuracy better than 
30 -50 % , the error is mainly due to a poor knowledge of multiparticle scattering 
FSI. Effects of soft FSI are as important as perturbative (and non-pertrurbative) 
QCD corrections up to the physical value of the b-quark mass, though they are 
parametrically suppressed in agreement with the color transparency (CT) picture 
of Bjorken |l]]. In a sense the situation is reminiscent of the case of the pion elec- 
tromagnetic form factor, which is dominated by the soft physics up to currently 
available momentum transfers Q 2 ~ 10 GeV 2 flT5"l|. There an additional suppression 



1/Q 2 of the soft contribution overcomes the leading hard contribution due to a large 
numerical factor. It may yet seem surprising that in our problem perturbative QCD 
is insufficient even at much higher momenta ~ m 2 B . 

A similar to our approach to the analysis of FSI effects in B-decay has been 
suggested in a recent paper || (see also p2| ). We differ from Refs. || |22| in two 



aspects. We concentrate on FSI reducible to the leading reggeon exchange, while 
the authors of || study worse understood elastic rescattering effects. Secondly, 
we calculate not only the phase but also corresponding corrections to the real parts 
of amplitudes. What if we apply the formulae of proceeding sections to the pomeron 
amplitude of elastic rescattering? We have checked that in this case we reproduce the 
conclusions of ||, i.e. the corresponding amplitude is not suppressed respectively 
to the direct one. Along with the two channel model used in || multiparticle 

1/2 

contributions do not seem to cancel the parametric m B dependence of the pomeron 
amplitude. It appears therefore that the pomeron breaks down the CT picture of Ref 
T[. We would like to mention, however, that assumptions underlying this conclusion 
are different from those done in the CT scenario, thus the final conclusion of Ref. M 
does not come too surprising. We feel that a space time picture of B-decays similar 
to those applied to high energy scattering would be desirable for a further study of 
this issue. 
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